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Abstract 

An energetic justification of a thermal component during inflation is given. 
The thermal component can act as a heat reservoir which induces thermal 
fluctuations on the inflaton field system. We showed previously that such 
thermal fluctuations could dominate quantum fluctuations in producing the 
initial seeds of density perturbations. A Langevin-like rate equation is derived 
from quantum field theory which describes the production of fluctuations in 
the inflaton field when acted upon by a simple modeled heat reservoir. In a 
certain limit this equation is shown to reduce to the standard Langevin equa- 
tion, which we used to construct "Warm Inflation" scenarios in previous work. 
A particle physics interpretation of our system-reservoir model is offered. 
PACS numbers: 98.80.Cq, 05.40.+j 

In press Phys. Rev. D 1996 

hep-th/ 960 11341 



1 



I. INTRODUCTION 



According to inflationary cosmology, the large scale structure of the present day uni- 
verse is essentially a kinematic outcome of exponential amplification of perturbing seeds in 
an initially smooth universe The differential microwave radiometer (DMR) on the 

Cosmic Background Observer (COBE) has made the first direct probe of the initial density 
perturbations through detection of the temperature anisotropics in the cosmic background 
radiation (CBR). These results are consistent with the scaling spectrum given by the in- 
flation model. They also reinforce previously known measurements, although done by less 
direct methods, that the amplitude of initial perturbations is 

A(k) = ^ = 10-=^ - 10-^ (1) 
P 

and approximately constant for all wavenumber k [^,0. Here p is the energy density in the 
present universe and Sp(k) is its rms-deviation at wavenumber k. 

We can understand the underlying kinematic origin of large-scale structure formation 
through a picture. Let us imagine observing the universe before inflation. It is a small 
patch that is growing at light speed. The process of inflation can be thought of as a rapid 
stretching of this patch in all directions. Comoving and physical coordinates are useful 
for further description. Comoving coordinates stretch with the patch so do not change for 
points that are stationary with respect to the patch. Physical coordinates express distances 
in terms of a physical measure, such as the local speed of light. A physical coordinate system 
is defined locally to a given point on the patch. It is useful sometimes to understand global 
distances in terms of physical units. For definitiveness, at the onset of inflation, let the 
comoving coordinates coincide with the physical coordinates. If we imagine the universe 
to be spherical with radius Rq at the onset of inflation, then during inflation, in comoving 
units, the radius remains the same whereas in physical units it grows as c^^Rq. Thus during 
inflation points on the preinflationary patch that are stationary in the comoving frame will 
rapidly move apart in terms of physical coordinates. 

Identify a point a on the preinflationary patch. Suppose that before inflation instruments 
are placed at a, which maintain communication with all points within reach by light since 
point a was created in the Big Bang. As inflation begins, communication is first lost with 
points furthest away from a. Immediately after inflation, only points that were initially 
very close to a will remain in communicative contact. To be definitive, just before inflation 
let the most distant points communicating with a be some distance Lq ~ 1/H. Here H is 
the Hubble constant with l/H = 5 x IQ-^^ GeV"^ = lO^^s cm = 3 x lO^^^ sec. The time 
interval is typically referred to as an e-fold. 

To understand the behavior of physical measures during inflation, let the origin of cosmic 
time, t=0, correspond to the beginning of inflation. The physical distance of a comoving 
interval Aa; at time t will then be e^^Axphys- On the other hand, for a light signal emitted at 
t = 0, computing along its geodesic ds'^ = 0, it will travel a physical distance (e"/if)(l — e"") 
after n e-folds of inflation. 

The implication of these two relations to physical correlations can be understood from the 
following example. For a signal emitted at t = 0, points fixed with respect to the comoving 
coordinates (comoving points) that are less than a comoving distance (l/if)(l — e~^) will 
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receive this signal before the first e-fold of inflation. However comoving points greater than 
a comoving distance 1/ H from a can never receive the light signal within the inflation stage. 
Only a sufficient time after inflation could such points communicate with a. In terms of 
physical distances, one concludes that comoving points that are greater than a physical 
distance at t = will lose communicative contact with a all during inflation. Since 
nothing is special about point a to any other point on the patch, on the large scale physical 
distances at t = that are greater than 1/ H will act incoherently all during inflation. Thus 
correlations larger than 1/H physical units are thereafter "frozen". This phenomenon is 
sometimes referred to as freeze-out or horizon crossing [Q. 

If a emits a second light signal at the end of the first e-fold, by the end of the second 
e-fold this signal again will have traveled a physical distance {e/H){l — e^^). However 
comoving coordinates will have stretched by a factor in terms of physical units. Thus 
only those points which are a comoving distance / H or less can receive the second signal 
at sometime within the inflation stage. Expressed in physical units, the conclusion stated 
above from the first e-fold repeats itself for the second e-fold: comoving points that are 
greater than a physical distance 1/H at t = 1/H hereafter will lose communicative contact 
with a all during the rest of inflation. Finally a signal emitted from a at the end of the 
n-th e-fold can be received by comoving points within a comoving distance e~'^/H from 
a. Thus comoving points greater than a physical distance 1/H at t = n/H thereafter lose 
communicative contact with a during the rest of inflation. 

Any local energy perturbations during inflation can affect a region of characteristic phys- 
ical length 1/H or less. The earlier a given perturbation occurs, the more elongated it will 
become due to inflation. As such the largest scales of energy density fluctuations in the 
post-inflationary universe arose from the earliest perturbations during inflation. 

In the standard scenario 0, inflationary expansion is assumed to occur within a super- 
cooled universe, in which the initial energy density perturbations were produced by quan- 
tum fluctuations. However energetics does not require a supercooled state. Furthermore, 
attempts based on this assumption have shown unnatural features. 

For inflation, naturality has played an important role. This is understandable since for 
phenomena that can not be directly observed, one attempts a description starting with the 
most natural expectations. The importance of naturality principles is to provide guidance 
from more familiar analogies with the hope of gaining predictability. For inflation we can 
understand naturality as both a macroscopic and microscopic one. Macroscopically we 
would like a description that rests with common-day experience. Microscopically it should 
be consistent with the standard model of particle physics. 

Under both categories, the standard scenario has shown unsatisfactory features [^0,[7[. 
Microscopically slow-roll scalar field dynamics requires an ultra-flat potential, although no 
such potential is required otherwise for particle phenomenology. Macroscopically, reheat- 
ing requires globally coherent radiation waves on the scale of the inflated universe. Local 
incoherent heat transfer is more familiar to experience. Furthermore, a globally coherent 
heating process requires a large scale radiator, which in the standard scenario is the inflaton. 
This raises the question about how the random inflaton field configuration before and during 
inflation attains quantum coherence at the end? 

We can accept that naturality principles for inflation do not have rigorous justifications 
and therefore can be abandoned or softened. However with this, any picture of inflation 
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based on familiar analogies would require modifications. This loosens theoretical constraints 
which otherwise are required to be consistent with only the limited data from observation. 
As such predictability from theory becomes less definitive. Nevertheless, if that is the way 
nature works, then that is the way it is. 

This would be an acceptable conclusion once all attempts for a natural explanation have 
been examined. If we digress back to this elementary point, we can rethink about the 
known ways to induce energy density perturbations. More general to quantum fluctuations 
are thermal fluctuations. We showed in that in the context of near equilibrium thermo- 
dynamics, during inflation thermal fluctuations could dominate quantum fluctuations under 
certain conditions. Neither energetics nor standard model dynamics precludes a thermal 
state during inflation. In [P] a model "Warm Inflation" scenario was considered for a self 
interacting scalar field. The solution had the desirable feature of satisfying observational 
constraints with minimal specifications placed on the field theory. The condition for slow- 
roll was shifted from requirements on the potential to a frictional force term that coupled 
the inflaton to a thermalized heat reservoir. By shifting to a frictionally produced slow-roll, 
it also gave a local heating mechanism. 

The first goal of this paper is to clarify the energetics in a thermal inflationary environ- 
ment. Secondly we derive from quantum field theory a generalized Langevin equation for 
the scalar inflaton field coupled to a modeled heat reservoir system. In certain limits, which 
we state, the equation reduces to the one we used in 0. Finally with these ingredients, 
an interpretation in terms of particle physics is attempted. To maintain physical clarity, 
our derivation is Hamiltonian based and performed in a cubic box with periodic boundary 
conditions. Appendix A provides a convenient reference to thermal properties of free fields 
that are used often in the text. It also relates our notation to standard form as well as 
shows how to take the infinite volume limit. Appendix B gives an alternative derivation of 
the Langevin equation from that in section (III). The purpose for this is discussed in section 
(III). 



II. ENERGETICS 

Let us account for the total energy in the inflationary universe. Consider a scalar inflaton 
field with Hamiltonian density defined with respect to the physical volume 

'H=\^' + V{<P). (2) 

First suppose that the inflaton is the only system in the universe. If during inflation the 
inflaton has negligible kinetic energy, (0^) ?s 0, then the energy density during inflation 
would be all potential from V^cpo). Such a situation could occur if 0o were at a local 
extrema where V^'(0o) = 0. Although the energy density remains constant in this case, the 
volume U of the universe in physical units would grow after n e-folds to U (n) ~ i^jje^" where 
U (0) ~ i?g is the initial volume of the universe just before inflation. In terms of total energy, 
Exin) = U{n) < 7i >, it would be after n e-folds of inflation 

ET{n) = V(0o)f/(O)e3" = ET{0)e^''. (3) 
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Turning to the more realistic situation in which there is kinetic energy, it is conceptu- 
ally helpful to first understand energy transfer in the classical limit for the inflaton field. 
The quantum mechanical problem is treated in section (III). For the classical limit in the 
expanding universe, the rate of change of the energy density can be expressed as 

^ = -3Hm (4) 

If in addition the inflaton field expels energy to some other system, one can express this, 
when treating the inflaton as an isolated system, by adding a dissipative term to the right- 
hand-side of eq. (^). If we choose the specific form to have the lowest time derivative and 
be even in the field, eq. (^) becomes 

dn{t) ^ _3^^2^^) _ f dt'^{t)V{t,t')^{t'). (5) 



dt 

This implies that the equation of motion for the inflaton is 

0(t) + m^{t) + 1* r(t, t')^{t')dt' + = 0. (6) 

In the next section we derive a quantum operator equation similar to this but also include 
a random force term and treat spatial variations. 

Having inserted a dissipative term in eqs. (^, energy balance implies that there must 
be some other system receiving this energy. If the second system is sufficiently large, it 
will act as a heat reservoir which induces fluctuations on the inflaton field. In the next 
section we examine a model heat reservoir system which we assume is thermalized. For the 
reservoir, we do not commit ourselves to a specific particle physics realization. However, 
in section (IV) a particle physics interpretation is offered. Furthermore, one may question 
the assumption of thermalization for the heat reservoir. More general would be some other 
statistical distribution. However for the present work we assume that the heat reservoir is 
in thermal equilibrium at some temperature T. Further treatment of this problem would 
require details about the dynamics beyond what we consider. Finally during inflation, the 
temperature also could be a function of cosmic time, but we consider it fixed. 

What is not an assumption and the important point to establish here is the energetic 
justification of the system - heat reservoir decomposition. The validity of this as well as 
consistency with the inflation solution requires 

(5p^(x,t) < p,(t) < p^(t). (7) 

Here 

5p<^(x,t) = 7^(0(x,t)) -7^(0o(t)) = \/'(0o(t))¥(x,t) (8) 

is the energy density contained in the fluctuations of the inflaton field 

50(x,t) = 0(x,t)-0o(t), (9) 

Prit) is the energy density in the heat reservoir, and p^if) is the vacuum energy density, 
with all of these evaluated at cosmic time t. 
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The first requirement, from tlie riglit inequality in eq. (|^), is tlie vacuum dominance 
condition needed for inflation and the second requirement, from the left inequality, is needed 
for the system-heat reservoir decomposition. In |^ the warm inflation scenarios that we 
found were for 

r = £l< 10-2. (10) 

This is sufficient to satisfy the first requirement of vacuum energy dominance. This state- 
ment can be strengthened for certain solution regimes found, which were consistent with 
observation and had R one to two orders of magnitude less than eq. (p!0|). The second 
requirement can be established from the observed amplitude in eq. (P since |]3|,p!0| 



lo-^-io-^^ r^ = « . (11) 



V P /he V0^+sP^ 



HC 



Here HC and HE mean horizon crossing and entry respectively, and k is the comoving 
wavenumber with horizon crossing occurring at cosmic time t such that |k|e-^* ~ H. On 
the left-hand-side 6p from eq. ([1|) has been equated to (5p0. In Q we found that the minimal 
thermodynamic requirement for a thermal scenario was ~ 0^, although for the model in 
1^ for all cases ^ 0^- In either case from eq. ([T^) we have 



10"^ - 10"^ (12) 



which satisfies the second requirement in eq. (|^). Thus energetically a thermal component 
can exist during inflation. 



III. FIELD THEORY 



In this section a field theory derivation of the operator equation of motion for the inflaton 
is given which has the form of a Langevin-like rate equation. In such a rate equation 
was postulated. Below we consider a simple model heat reservoir system which is coupled 
linearly to the completely interacting inflaton field. The total system-reservoir Lagrangian 
is 



Lt — Ls + + L/, 



(13) 



where on the right-hand-side the Lagrangians are Ls for the inflaton (system), L/j for the 
reservoir, and Lj for the interaction between system and reservoir. The inflaton's La- 
grangian, Ls, can have an arbitrary potential and accounts for the expansion term. It 
has the familiar form 



3Ht 



2((^o0(x,t))^ 



(e-^*V0(x, t)f - m'<p\^, t)) - V(0(x, t)) 



(14) 



where the potential can have the general expansion 
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n0) = f:§0"(x,t). (15) 

n=3 

In this paper we derive the effective operator equation of motion for the inffaton, but we do 
not study it any further. Thus we will not address the issues of renormalization, which would 
be connected with solving this equation. As such, our formal derivation is valid whether Ls 
represents an elementary or effective field theory. 

The heat reservoir Lagrangian, Lpt, is modeled as a set of free fields each characterized 
by a mass fii. It is written as 



i 



E //'xe^^*^ [i^oV^Mr - (e-^Vr^.(x,t))^ - /ik'(x,t)] . (16) 



Each reservoir field is coupled linearly to the inffaton with a coupling constant through 
the interaction Lagrangian 

^/ = - E «^ / c?'xe3^*r7,(x, t)0(x, t). (17) 
. Jv 

Note that the coupling constants {aj} carry engineering dimension 2. 

The derivation given below treats the energy transfer between the system and reservoir. 
However the modeled reservoir in eq. (|16D does not have internal interactions. Furthermore 
energy transfer between the inffaton and the reservoir is not fully treated. To understand the 
latter two points, note that our problem has a difference to standard problems in Langevin 
dynamics. The reservoir acts as a large system on the fluctuations of the inflaton field, 
which is a standard situation for applying Langevin dynamics. However the difference 
is that the vacuum energy of the inflaton, and particularly through the zero mode, acts 
as a energy source for maintaining the energy density of the reservoir, which otherwise 
would diminish due to inflationary expansion. This introduces two complications. The flrst 
concerns thermalization for the reservoir. In a standard Langevin problem, it is assumed 
that the system will weakly interact with the subsystems of the reservoir. There can be many 
subsystems so that the total effect of the reservoir on the system can be strong. However 
each subsystem of the reservoir is affected only weakly by the system. In such a circumstance 
the issue of maintaining an initially thermalized state for the reservoir is not acute. As such, 
internal reservoir interactions are not crucial to know. In our problem we can accept that 
the couplings in eq. (|T^ are small (although the derivation below holds for arbitrary 
aj's). The question of concern is for whatever vacuum energy that is transferred into the 
reservoir, can it thermalize on a time scale shorter than 1 / H7 Even if our modeled reservoir 
had internal interactions, to prove thermalization from flrst principle dynamics would be 
complicated (some attempts in one-dimensional models are given in ||TT[|). 



In our treatment we assume that thermalization occurs and make appropriate by-hand 
adjustments. These are detailed in the derivation, when relevant. In the next section we will 
return to the elementary question of thermalization again. Alongside with this problem is 
the second complication, which is in regards to energy transfer from the inflaton vacuum to 
the reservoir. Our derivation accounts for only sub-horizon scale physics. As such it treats 
the fluctuations induced on the modes of the inflaton fleld by the reservoir, while the modes 
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are sub-horizon scale, kphysicai > H. However our derivation is not justified for treating tlie 
zero mode interaction between tlie inflaton and tlie reservoir, since tliat also involves super- 
horizon scale physics. The suggestive guess is that the "naive" operator equation derived 
below for the zero mode, but without a thermalization assumption on the super-horizon 
scale modes of the reservoir, is approximately valid. We will return to this issue in the 
conclusion. 

It is worth noting that in certain limits the thermal inflation problem reduces to the 
standard Langevin problem, although this limit did not prove useful in the warm inflation 
scenarios we considered in [0. This limit is given in Appendix B as well as the derivation 
of the Langevin equation in this limit. 

Our first goal is to derive the effective equation of motion for the inflaton with the 
reservoir field variables eliminated. The equation is valid for a time interval |At| < H~^. 
We will derive the equation for an arbitrary comoving mode of the inflaton field and for 
the m-th time interval, tm to tm+i, with m arbitrary, where tm = mAt and to = is fixed 
as the starting time of inflation. A complete solution for the inflaton's evolution can be 
obtained by piecewise construction over all time intervals. The formal derivation is not 
different for physical wave- numbers that are sub-horizon or super-horizon scale, although as 
stated above, the approximations leading to the derivation are valid only for the former. 

In solving the equations of motion for the reservoir fields, three approximations are made 
in the m-th time interval for every m. Firstly the red-shifting factor of a given comoving 
wavenumber k„ is held fixed at e"^^*™. Secondly, at the beginning of every time interval, 
the state of the reservoir field operators is readjusted. Finally the uncoupled modes of the 
reservoir fields are assumed to obey a canonical distribution with respect to their physical 
frequency 

Here time t has been demoted into the subscript to signify that it is treated as an adiabatic 
parameter as far as the wavenumber is concerned. 

The first approximation is made to simplify the calculation so that it can be solved ana- 
lytically. It can be dropped if one is willing to apply more sophisticated methods of solution. 
The latter two approximations are physically motivated. They are by-hand treatments of 
the interactions amongst the reservoir fields. Details about the second approximation are 
given within the derivation when relevant. The third approximation implements our ther- 
malization assumption for the reservoir. In the derivation we are careful to separate these 
two approximations. The second is made at the operator level. The third is a statement 
about the state. For it, we are firstly assuming that the description of the reservoir is statis- 
tical and secondly that the particular distribution is canonical. In the next section we argue 
that the reservoir state is created from quantum decay processes, in which case a statistical 
description is inherently required. The assumption of being a canonical distribution seems 
the most obvious first guess. 

The final equation of motion for the inflaton will be stochastic since the reservoir state 
is specified by a statistical distribution. The final equation of motion for the inflaton su- 
perficially will appear nonconservative since the reservoir fields are going to be eliminated. 
Under certain conditions placed on the reservoir Hamiltonian, we derive the Langevin equa- 
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tion used in M. This limit is examined and we verify the fluctuation-dissipation theorem in 



its standard form 12 



The derivation below follows well known methods from non-equilibrium statistical me- 



chanics which have been refined over the years [Il3|-[T6[]. A primary motivation for this 
ongoing effort has been to understand the universal properties of the Langevin equation and 
to obtain a possible explanation from first principles [p!3| , p!7| . This believed universality is 
one reason for us to start our study of stochastic dynamics for thermal scenarios with the 
Langevin equation. Our derivation below follows closest to |jl3[ and the model for the heat 



reservoir follows [|T^]. We have made some modifications to these works, which were for 
quantum mechanical models, in order that we can treat a quantum field and account for the 
expansion term. In Appendix B there is an alternative derivation which obtains the rate 
equation in certain limits for an arbitrary reservoir Hamiltonian. 

Before proceeding. Let us review literature that is related to the present work. Applica- 
tions in cosmology using Langevin dynamics have been done in although both our 
methods and motivation differ from these works. Studies of finite temperature field theory 
in Robert son- Walker Universes have been done in |2^. Finally a calculation with similar 
objectives to ours in this paper is given in [^]. There a path integral derivation is presented 



of the infiaton's evolution when coupled to a thermal bath in de Sitter space. However the 
authors did not completely examine the dissipative properties of such a system. As such 
they apparently missed the connection to warm infiation type scenarios, which for us is the 
starting motivation to the present formal exercise. 

We perform our derivation in a cube for the three spatial directions, which is centered 
at the origin with sides at ±L/2 in each direction. The Fourier expansion of a generic field 
is from eqs. (|A9|) 

x(x,t) = ^Ex(kn,t)e^^"■^ (19) 
where throughout this paper we use the notation k„ = 27rn/L, with n = {nx,ny,nz) and 

oo oo oo 

E E E • (20) 



" na;=— oo 71^=— oo nz=— oo 



The argument of the coordinate (x) and the momentum (k) space fields are always given to 
distinguish the two. 



The system-reservoir Lagrangian in terms of the Fourier modes as defined in eq. (|T9|) is 
Lt = Ls + Lji + Lj, where 



I E [<^(kn, t)0(-k., t) - (e-2^*k2 + m')<j>{K, t)0(-k„, t)] - Vp] , (21) 



^3Ht 



L3 



E E ^ [viiK, t)rji{-K, t) - (e 2^*k2 + ^xl)r]i{K. t)r?^(-k„, t) 



(22) 



and 
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L3 



with 



^3Ht 



7rx(kn,t) 



-xi-K,t). 



Converting to the Hamihonian Ht = ir^^' + J2i T^rnVi ~ -^t, we obtain 



,3Ht 



(23) 



Vp = L' d^^V{4>) = E 7^ E ■ ■ ■</.(k^„_„t)</.(-k„, k^_,,t). 

n=3 mi,--,m„_i 

(24) 

The conjugate momentum to any scalar field x(k„,t) is 

dL 



e-^'^'L^iK, t)n^i-K, t) + -j^ie'^'^'kl + m')0(k„, t)0(-k„, t) 



(25) 



(26) 



^^ = EE^ 



(27) 



and Hi = —Lj. Our notation is that all Hamiltonians have some specifying subscript, which 
then leaves the Hubble constant to be H. 

To quantize the theory the postulated equal time commutation relations (CCR) are given 
in eqs. ( [A^ ) and ( |A11| ). The operator equations of motion from Ht are 



(28) 



TT^iK, t) = i[HT, Tl^iK, t)] 



^3Ht 



L3 



e-2^*k2 +m^)0(k„,t) + 



6Vf 



5(f){-kn,t) 



+ E"i^i(kn,^) 
i 

(29) 



ri^iK,t) = e-^'''L\X^n,t) 



and 



.,3Ht 



L3 



e ^^*k^ + ^il)r]i{kn, t) + a,0(k„, t) 



(30) 



(31) 



The resulting second order field equations are 



(k„, t) + 3//0(k„, t) + [e-^^*k^ + m^]0(k„, t) + 



5(j){-kn,t) 



+ Eai^i(kn,t) = (32) 



10 



and 



UK, t) + 3HT]i{K, t) + [e-'^*k2 + /i,"]r/i(k„, t) + ai0(k„, t) = 0. 



(33) 



We now implement our second approximation discussed earlier. At the beginning of 
every time interval 1^-, we readjust the state of the reservoir fields as 



m— 1 J 



r/-(k„,0;t. 



(34) 



where ?7i(k„, t; tm-i) and //f (k„, t; tm) are the solutions of eq. (0) at time t respectively with 
and without the coupling term to and with the frequency ^n{t) held fixed at respectively 
tm-i and tm- Formally this operation can be viewed as a set of impulsive forces that act on 
the reservoir. The first purpose of these adjustments is to add sufficient energy so that the 
reservoir's energy density with respect to physical volume remains constant. Secondly these 
adjustments are an external treatment of interactions within the reservoir. They shift each 
field back to its free field state at t = but with its physical frequency decreased by a little. 
This operation is in preparation for the thermalization assumption we will make below on 
the reservoir's state. The first purpose given above acts as a constraint on the infiaton's 
evolution. Consistency with the infiaton's equation of motion for a particular model was 
demonstrated in |p. The discretized treatment simplifies the calculation. One expects the 
actual dynamics to be smooth and continuous. 

The solution for the oscillator field rji for the time interval [tm,tm+i] is 



Qi2 



-3H(t-tm) 

(k„,t) - 0(k„,tm)e 2 cos^2^, - 



-e-t^(*-*-)/ rft'cosfi^(,„)(t-t')e^ 



^H{t'-t^) 



3H 

[k^,t') + —<p{K,t'] 



(35) 



where 



-2HtV2 



9H^ 



(36) 



and 

Vi (kn; t] tm) = 6 



3Ht 
2 



vKK,o) cosni^,^^{t) + 



— f^r^^{k^,0) + L\XK,0) ] smQ\,.{t) 



(37) 



is the solution for the free reservoir fields with the time dependent physical frequency eq. 
(p!8|) of the comoving mode held fixed at tm- As in eq. (p!8[) , the subscript n(t) in eq. (|36D 
refers to the physical wavenumber of the comoving mode k„ at cosmic time t. Substituting 
eq. (|35|) into eq. (0), we obtain for tm < t < tm+i 



{K,t)+3H<p{K,t) + e~ 



\H{t-t^) 



E 



a,; 



dt'ei^(*'-*-)cosl]t,(,^)(t-t')0(k„,t') 
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+ (e-^^*k^ + m2)0(k„,t)-^ 



ai 



n 



n{tm) 



0(k„, t) - <f)(kn, t„)e-i^(*-*'") COS {t - tm) 



^-lH{t-tm) 



+ 



3^ 
50(-k„,t) 



where 



(38) 



(39) 



Up to this point, no statistical assumption has been made. Following our earlier discus- 
sion, we now assume that the free reservoir fields are canonically distributed. The statistical 
mechanics for the reservoir system based on our above approximations is the same within 
each e-fold as in fiat-space. For this Appendix A has been provided as a useful reference. 
As a clarification, no approximation has been made in treating the infiaton system. For the 
reservoir fields from eq. ( |A19| ) for tm <t,t' < tm+i, we have 



{{ViO^njt tmi ^m) ViO^n^^ tm]tm))T 



_p-^(i+i'-2tm)x 
2 ^ On-n' 



+ 



+ 



1 



90* 2 



A, ,t ■ n(tm) 



COSCiJ, 



.(t™)(^-0-cos^;(,^)(t + t'-2t^) 



3/7 



sin<t„)(i + ^'-2t. 



coth 



n(tm) 

2T 



n{tm) 



(40) 



Eq. ( |38D with the assumption eq. ( ^OD is the general rate equation and the main result of 
this section. For use later, the symmetric correlation function of the force operator is 

ST{K,t,t';tm) = ^{{r]{kn,t;tm)r]i-K,t';tm) + r]{-kn,t';tm)v(kn,t;tm)))T 



on. 



n{tm,) 



+ 



+ 



Qi 2 



cos^;(i^)(t - t') + cos<(,^)(t + t' - 21, 

- t') - COS^;((^)(t + t' - 2tm) 



3H 



sinc^;(t„)(t + t'-2t^) 



coth 



n{tm) 

2T 



(41) 
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We next derive a limiting form of eq. (pSl). Consider a large number of oscillator fields 
rji that are represented by a continuous distribution as 

E ^ / df^Nif^), (42) 

i 

where N{fi) is the spectral weight function. Assume that N{fi) is nonzero only in the 
interval fii < fi < fi^ with = fi^ — fii and the corresponding definition A^^^^^ = 
e-2mk2 + ^2 _ Mi _ ^6-2^^2 +^2 _ 9|i^ ^jg^ ^ssumc that 

= f (43) 

n{tm) 

is a constant over the interval A^. If A^^^^^ ^ H, then for l/A^^^^-, <^t<l/H 

r d^N{^){^fcosn^^^,J = 2T6y^. (t). (44) 

The subscript on the right-hand-side is to indicate that the 5-function is smeared over a 
time interval l/A^^^^-j. Inserting eq. ( PD into eq. ( |55D we obtain 



0(k„, t) + (SH + r)0(k„, t) + (e-2^*k^ + ^ + ^)0(k„, t) 

TT Z 



+ 1 + Sre-^^^*-*™ V(kn, tm)5uA» (t - tm) = v'iK, t; trn) (45) 
()0(-k„,t) ' "Cm) 

where 

V(k„, t; tm) = - j dfiN{iJ,)a{fi)r]^(k„., t; tm] yu). (46) 

In the limit of zero expansion, = 0, we obtain from eq. (l45|) 

2rA 

0(k„, t) + r0(k„, t) + (k^ + ^)0(k„, t) 

TT 

+ . , . / + 2r0(kn,t^)(5i/^n (t - 1„) = r7o(k„,t) (47) 

where 77Q(k„,t) is as in eq. (^) but with if = in eq. (^). Aside from the last term 
on the left-hand-side, which is nonzero only for a very short time interval, 1/A^, this is the 
standard form of the Langevin equation and the type we used in 0. In the limit T ^ oo 
we find from eq. (^) and eq. ( PD 



^((r/o(k„, t)r/o(-k„, t') + r7o(-k„, t')^7o(k„, t)))T^oo ^ 2L=^rT5(t - t'), (48) 



which verifies the fluctuation-dissipation theorem ||12 
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IV. INTERPRETATION 



In the last section we obtained the generahzed Langevin equation (^Sf), which induces 
fluctuations on the inflaton as well as frictionally damps its motion. To obtain this equation, 
we had to introduce a heat reservoir of light particles that interacted with the inflaton. In 
P one form of rate equation (0), namely eq. {^7\), was shown to be numerically successful 
for inflaton dynamics. This evidence provides sufficient motivation to seek an interpretation 
of the heat reservoir in the context of particle physics. The heat reservoir of the last section 
is a toy model. Here we want to think generally about a second system of light particles 
which acts random force on the inflaton. 

For the dynamics described in the last section to be realizable in the real world, in par- 
ticular two properties are needed for the reservoir. Firstly there must be some mechanism 
available during inflation to produce the reservoir particles. Secondly these particles must 
interact rapidly on the scale of the expansion time 1/H. In section (III) we assumed that 
they could thermalize. This stringent of a condition is not needed for eq. ( pS] ) to describe 
stochastic evolution of the inflaton. Nevertheless, we keep with our thermalization assump- 
tion and will see how close we can come. In addition to these two properties, a less important 
point of interest is what direct interpretation the explicit reservoir Hamiltonian (0) in our 
simple model has. 

To address these issues, we first make three general observations. Firstly, the zero mode 
of the inflaton is evolving quasi-statically during its roll down the potential. Secondly, 
quantum processes will occur within micro-physical scales, so for physical distances and 
cosmic time less than 1/H at any e-fold of inflation. The third observation is based on 
our analysis in 0. There we found that the energy density of the heat reservoir was close 
to the GUT scale Mq^t- particular the temperature range of the thermal scenarios we 
considered that were consistent with observation was T ~ (0.01 — 0.03)Mgut- This implies 
for a single light species a corresponding energy density in the range 

Pr ~ (10-« - 10-6)M^uT- (49) 

For Mgut = 10^^ GeV= 1.782 x 10^^ g and for the lower bound on this means 

Pr > 10-^M^uT = 10^2 GeV^ = lO'^^g/cm^ 

~ lO^^'^Pncutron star ~ lO^'^PHI-coUisionj (50) 

where Pneutron star and Pffl-coUision are the energy densities for respectively a neutron star and 
an upper bound estimate for a heavy-ion collision. For our present purpose, we make the 
general observation that the energy density of the heat reservoir is large on the scale of the 
largest probed energy densities. 

Observation one provides suitable conditions for particle creation from quantum decay 
of the inflaton similar to pair-production in a strong quasi-static electric field. This gives 
a mechanism for producing the particles in the heat reservoir. An outcome is that the 
description of the reservoir's state must be inherently statistical. 

Having made a particle association for the heat reservoir, for completeness of our model 
in section (HI), we make an interpretation of Hn in eq. (pTj). The most direct statement is 
that Hji, with an appropriate spectral weight function eq. (^2]), is a general representation 
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of scalar quasiparticle excitations in the noninteracting limit. On the other hand, thinking 
in terms of elementary fields, if we apply dimensional counting to Hj in eq. (P) we find that 
the reservoir fields rii{'x,t) carry engineering dimension 3. This suggests that in our simple 
model each reservoir field represents a fermion-antifermion pair (ie. qq, ee etc.), which 
interacts with the inflaton through a Yukawa coupling. Hji then has the interpretation of 
an effective free field Hamiltonian for paired // states. In QCD these would be the simplest 
color singlet states constructible from quarks. Beyond these "kinematic" statements, we can 
not say more at present. Almost nothing is known about matter at the energy densities of eq. 
(0). We therefore must leave unanswered to what extent Hr in eq. ( |77| ) and the spectral 
weight function in eq. ( ^31) approximate a heat reservoir system which can be derived from 
a realistic particle physics Lagrangian. 

Returning from this digression, next we will examine what expectation there is for ther- 
malization in the reservoir. For this let us review the situation up to now. Eq. (^) contains 
the dynamics that governs energy transfer between the reservoir and the inflaton. The full 
internal dynamics of the inflaton system in isolation are also treated. Thus eq. ([38| ) is valid 
for an arbitrary scalar inflaton potential and for reservoir fields that do not interact directly 
with each other. We discussed in section (III) that in the rapid expansion environment 
during inflation the internal interactions of an actual reservoir system are essential for dis- 
tributing any influx of energy. Having established above that the state of the reservoir is 
inherently statistical, knowledge of internal reservoir dynamics is in particular important for 
addressing the question of thermalization. In the last section we got around the problem by 
assuming that thermalization occurs and then made suitable adjustments to our equations. 

To check the assumption of thermalization is a difficult problem even if one is given more 
dynamical details about the reservoir. However based on observations two and three we can 
get oriented with the scales involved in the problem. For this let us consider an unjustifiable 
classical model of light hard rocks of average energy E, center-of-mass two-body elastic cross 
section a, and at a energy density procks of order eq. (^Dp. The mean free path in this system 
is 

E 

Crocks ; (^-^) 

Procks'^ 

which for relativistic rocks implies a collision frequency t^rocks = l/4ocks- 

We assume that the elastic cross section is independent of energy. Up to logarithmic 
corrections, this corresponds to Regge asymptotic behavior from single pomeron exchange 
with slope one (for a review please see |^ and references therein) . For a conservative lower 



bound to elastic pp-scattering at high energy, we use aei = 10 GeV^^. The energy of a 
typical light particle at temperature T is ~ T. Taking the lower bound from eq. (|50D 
Procks = 10^^ GeV^ which corresponds to T = O.OIMgut, we find 

t^rocks = 10^° GeV ^ 10^° i/. (52) 

It is known that most of the elastic cross section a^i is increasingly in the forward direction 
with increasing energy. Accounting for this, we take the average momentum transfer to 
be \t\ = 0.025GeVl Let us assume statistical independence between collisions for a typical 
light particle. Then within one e-fold, the root-mean-square energy that the particle transfers 



15 



will be about an order of magnitude greater than its total energy and it samples all angular 
directions in phase space. Optimistically the range for eq. (0) could be increased by a 
factor 10^ by using the upper bound on Procks in eq. ( ^Q]) and increasing cXez by over a factor 
10. This range presents a case that some sort of GUT-plasma could thermalize within an 
e-fold at GUT-scaled energy densities. 



In section (III) we derived a quantum Langevin equation, but we can think of it as a 
classical Langevin equation. For this we can make a picture of the dynamics for a classical 
infiaton field 0'^(x, t). We define the order parameter in the classical system as 



and it will sometimes be denoted as (j)^. For the classical theory we will make a mechanical 
model and then associate it with the classical field. We will work in the infinite volume limit 
where correspondence between the classical and quantum fields can be stated more simply. 

Let us first recall some details about the infiaton Hamiltonian. In thermal scenarios the 
symmetry breaking potential V{(j)) can be of the standard double-well type without an ultra- 
fiat region as demonstrated in 0. The relevant magnitudes which define the potential are the 
same as in the standard scenario. Thus the vacuum energy V{Q) ~ Mgyrp, the minima of the 
wells are at ~ i^fcuT, and the curvature at the minima are V^"((^m) ~ ^gut- Inflation 
begins when the infiaton tunnels out of the meta-stable minima at = to < 0i < (^m- 
We will define this moment as the origin of cosmic time t. Our classical analogy begins at 



With these preliminary remarks, let us state our classical analog model. We picture the 
classical scalar field as composed of heavy particles of mass = \JV"{(f)m)- They sit in a 
cubic lattice arrangement with massless (anharmonic if preferred) springs connecting them. 
Kinetic energy is associated with the particle masses and their motional energy. Potential 
energy is associated with the energy in the massless springs. The order parameter (j)'^{t) is 
interpreted to be related to the number density n<^(t) of 0-particles. To obtain a precise 
relationship we will equate the order parameter with the vacuum expectation value of the 
quantum Heisenberg scalar field at t=0 as 



We will assume that the vacuum is a coherent state of zero momentum particles. Using the 
number density operator 



V. MECHANICAL ANALOGUE 




(53) 



t=0. 




(54) 




- hm V a. 

' n 



(55) 



we arrive at 
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In eq. (^5]) () — lim means evaluate the operator matrix element in finite volume and then 
take the limit V —>■ oo. In our classical model we define the order parameter at all time with 
this relation, 



If one prefers, this can be regarded as an unmotivated definition, which relates the the 
classical mechanical model with our classical field model. 

The initial conditions at t = are the order parameter is at (f)'^{0) = (pi, the masses 
are at rest with number density given by eq. (|57D , and the springs are elongated such that 
the cubic symmetry is maintained but the fundamental cell is bigger. When the system is 
released a breathing mode commences, in which the fundamental cell contracts to a minima 
and the expands back to its original size at t = and so on. Correspondingly, the number 
density oscillates. 

To this system, we will add many light ?7-particles of mass -C that are thermalized 
at temperature T. This system forms the heat reservoir. The r^-particles randomly will 
impart momentum on the 0-particles. This will perturb the breathing mode of the 0-particles 
as well as excite transverse modes on the 0-lattice. The system so described is similar to 
that for Brownian motion, except now there are several Brownian particles and they also 
have an interaction amongst themselves. With minor changes, one can understand the rest 
of the statistical mechanics for our system from Brownian motion. Slow-roll motion in our 
classical model corresponds to over-damping, in which the fundamental cell slowly contracts 
to equilibrium with no oscillations. In the process the potential energy in the 0-system is 
transferred to the heat reservoir. One can also expand this picture to include infiationary 
expansion and ?7-particle creation, but we will leave that for the reader's thoughts. 



VI. CONCLUSION 

In section (III) we obtained the rate eq. ( PS]) which describes the dynamics of the 
infiaton's comoving modes while they are sub-horizon size. For the zero mode, eq. ( pH]) is 
therefore unjustified, since it does not treat super-horizon scale physics. What one needs 
is an equation that describes the transfer of potential energy from the infiaton's vacuum 
during its roll down to the potential minima. In our model, the energy is transfered not 
only into the infiaton's kinetic energy, but also into the fields of the reservoir. 

In a local patch and for a short time interval, eqs. (^) would be valid for all modes 
which are ad-hocally defined with respect to the patch. A network of such patches could 
be constructed to cover the infiating universe for some chosen cosmic time interval. In 
each patch the respective Langevin equation could be solved with solutions matched at 
the boundaries. To compute the initial sub-horizon scale energy density perturbations, 
recall that they are determined by the sub-horizon scale modes. Therefore this computation 
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could be done from the solution of the Langevin equation from any one patch. For the 
zero mode, there are two points of interest. First we want to know the global energy loss 
by the infiaton. Due to the lack of coherence between patches, the expectation value of 
the infiaton (order parameter) need not be the same in different patches. Different local 
regions of the universe could then be in different stages of inflation, depending on the local 
order parameter. For computing the total change in vacuum energy, one could sum the 
result from the independent patches. However turning to the second point, for inflation 
to solve the smoothness problem, one requires the order parameter to be approximately 
globally homogeneous. One may argue that the local conditions in all patches would be 
similar, so that the order parameter should not differ by too much in different patches. 
To verify this would require further investigation within specific models. In particular, to 
meet observational constraints, one wants to know the degree of homogeneity of the order 
parameter within an inflating region for a duration of at least 70 e-folds. All statements so 
far apply equally well to scenarios based on both quantum and thermal fluctuations. There 
is one difference. Thermal scenarios, do not require the quantum state of the inflaton's zero 
mode to be coherent over the entire inflationary patch. This is required in the standard 
scenario for reheating. 

The derivation in section (111) is formally valid for a fundamental or effective infiaton 
Lagrangian Ls in eq. (0). Given an initial fundamental Lagrangian for the infiaton, one 
anticipates that to formulate a consistent stochastic process, it will require high frequency 
interactions of the infiaton with itself and with other fields to be treated first and appear 
in the form of an effective Lagrangian Ls in eq. (0). We have not developed in this paper 
any methods for obtaining the appropriate effective Lagrangian one should use for Ls- 

In we assumed near equilibrium dynamics applied everywhere. We therefore used 
the finite temperature effective potential in Ls- However nonequilibrium considerations 
such as done in ||2^ for flat space may be important to extend for the inflationary era 
of expanding space. Generally for a scalar infiaton field coupled to gauge fields, such as 
in GUT theories, one should expect quantum corrections from self-coupling and quantum 
plus thermal corrections from coupling to the gauge fields to give a modified potential in 
Ls- Since the relevant temperature scale for slow-roll dynamics is of order or less than the 
infiaton mass, one does not expect thermal corrections from the inflaton's self-coupling to 
be substantial. 

To summarize, in this paper we showed that the presence of a sizable thermal component 
during inflation is energetically allowed. Following our hypothesis in [Q, we treated the 
thermal component as a heat reservoir, which interacted with the infiaton in the manner of 
Langevin dynamics. We demonstrated in section (HI) that the Langevin equation postulated 
in p could be derived from quantum field theory. In we showed that such a dynamical 
arrangement, for the specific case of a double-well scalar potential, was numerically successful 
in describing observational data without fine tuning the coupling constant. 

One of the improvements suggested by our calculations there was to have a time de- 
pendent temperature that smoothly varied by a factor 2-3 over the course of the roll-down 
period. Not only was this numerically preferable, but it seems that thermodynamics does 
not restrict some variation in the temperature just before to just after inflation. As such 
the general dynamical equations during inflation should be able to accommodate such a 
variation, although specific models may make specific predictions. 
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In the reservoir was defined by its thermodynamic properties, but was not given 
a dynamical representation. Thus we had no underlying mechanism that could explain 
temperature variations. In section (III) our model supplies an explanation. However it is 
clear from our treatment there, that it was an arbitrary decision to hold the reservoir's 
energy density constant with respect to physical volume all throughout inflation. Just as 
well within the same dynamical framework, it could increase or decrease and in different 
time periods do either. However this is as far as the model in its present form can go. On 
the plus side, the model goes far enough to dynamically explain a time dependent reservoir 
energy density. However the model is not sufficiently constrained to do better than this. The 
options to improve this situation are to determine a further constraint from either formal 
means or phenomenological motivation. 
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APPENDIX A: STATISTICAL MECHANICS OF FREE FIELDS 

Some statistical mechanical properties of free field, which are often used in the text, are 
given here. We consider a free field Hamiltonian defined in a box of volume Vl centered at 
the origin with sides at ±L/2 in all three spatial directions. We give results in the discrete 
form in Vl for finite but large L and then in the infinite volume limit, V^o, where L — > oo. 
The free field Hamiltonian which we study is 

H = f d^^l [vr^ (x, t) + (Vx(x, t)f + m\\^, t)] + C, (Al) 

where C is a constant that later will be set so that the lowest energy state is at zero. The 
fields satisfy the equal time canonical commutation relations (CCR) 

[X(x,t),x(x',t)] =0 
[7rjx,t),7r^(x',t)] =0 

K(x,t),x(x',t)] = -z5(3)(x-x'). (A2) 
The field have the expansion in terms of creation and annihilation operators 

1 °° 1 

X(x, t) = —j^ -=[a„e^(''"— *) + at e-^C'"— *)] 

^ n,,ny,n,=~oo V ^^n 

OO j 

vr,(x, t) = ^ Ke^^'^"-™*) - 4 e-^(^— *)] (A3) 

where 
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27rn 



(A4) 
(A5) 



UJ. 



and from eq. 



[a„,a„/] = [a],,ajj,] = 

[«n, = 5nn'- 



We define the Fourier transform of the fields as 



7rx(k„,t) = yr / (i^x7r(x,t)e 



ikn-X 



(A6) 

(AT) 
(A8) 

(A9) 
(AlO) 



From eq. ( [A^ ) one can verify that the Fourier space modes of the fields satisfy the equal 
time CCR 

[x{K,t),x{K',t)] = o 

[7r^(k„,t),7r^(k„/,t)] = 

[7r^(k„,t),x(-k„',t)] = -iSnn'- (All) 
The expansion of the Fourier-modes in terms of the creation-annihilation operators is 

r3/2 

xiK,t) = ^=[a„e-*^"* + al,e*-"*l 



The constant C in eq. (|A1|) is be defined as 

1 



C 



(A12) 



(A13) 



Two equivalent forms of the Hamiltonian H in eq. (^T|) are given below. Writing 

H = Y.Hn = Y.Hk„^ (A14) 



where in terms of the creation-annihilation operators in eq. (|A8| ) 
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Hn = UJnal^an 

and in terms of the Fourier- modes in eq. ( A12 ) 

Hk„ = Y%(kn, t)7T^{~K, t) + — (k^ + m^)x(k„, t)x{-K, t). 



(A15) 



(A16) 



The partition function, Z = Ti {e -^/^ j , evaluated in the creation-annihilation basis is 
a product Z = Yin ^n, where 



Z„ = T.(e-«../-) = ^-^. 



(A17) 



The two-point thermal correlation function for the creation-annihilation operators is 
^Tr (e^4a„/) = ((4a„/))T = Snn' -I _ 



((anan'))T = ((44'))t = 

and for the Fourier-modes is 
^Tr(eT^x(kn,t)x(k;,t')) = ((x(k„, t)x(k;, t'))T 



(A18) 



((7r^(k„,t)7r^(k;,t')))T = 5n,^ 



'21? 



coth 


( 




K2t) 


coth 






K2t) 



coth ( sinu;„(t — t') + icosuJn{t — t') 



From the last relation note that 

(([7r;^(k„,t),x(k^,t')]))r = -z5„_„/ cosu;„(t - t'). 
The total ensemble averaged energy of the thermalized system is 

UvAT) = T.U^ST) 

n 

where 



(A19) 



(A20) 



(A21) 



(A22) 



Let us examine the infinite volume limit. The creation and annihilation operators must 
be rescaled as 



c(k„,) = 1 77- ) a„ 



^.3/2 

27r/ 



(A23) 



and similarly for aj^. Substituting eq. (|A23[ ) into eq. ( |A3| ) and identifying 
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L 



L^oo 



L 



we get the continuum form 

[(27r)32l]V2[<k)e 



+ cHk)e 



'i(k„-a;-c<;k)*)l 



(A24) 



(A25) 



and similarly for -K^ip^^t). For the momentum space Fourier transforms of the fields in the 
infinite volume limit, Vl must be replaced by Vqo in eqs. ( |A9|) and ( |A1CI| ) and 7r^(k, t) must 
be defined without the factor l/Li^. The Hamiltonian in the forms eq. ( |A15| ) and eq. (|A16|) 
in Voo are respectively 



H= / rf=^kcukc^(k)c(k) 
rf^k 1 r 



(27r)3 2 

The energy density in Voo is obtained as 



7r^(k, t)7r^(-k, t) + (k^ + m')x(k, t)x(-k, t) 



From eq. ( [MID and (^22|) we find 

«(T): 



(27r)3e'^k/T 



(A26) 



(A27) 



(A28) 



The familiar formula for the energy density of electromagnetic radiation is twice this, due 
to two polarization states for the photon, and with m = so that cj^ = |k|. 



APPENDIX B: ALTERNATIVE DERIVATION OF THE RATE EQUATION 

For a system with a nearly harmonic potential, such as within one well of a double-well 
potential, one can derive the Langevin equation with a general reservoir Hamiltonian Hr in 
the regime 

< r < u;, (Bl) 

where H is the Hubble constant, F is the dissipative constant and u) is the frequency of 
the system with no anharmonic terms in the potential. We perform the derivation in this 
appendix. It follows the method of [|14|. As in section HI, it is done in a cube for the three 
spatial directions, which is centered at the origin with sides at ±L/2 in each direction. 
The Hamiltonian we examine is 

Ht = Hs + Hr + Hj, (B2) 

where Hs is the system Hamiltonian given for the inflaton field in eq. (]26|), if/j is the 
Hamiltonian for the reservoir system when it is uncoupled to the scalar field system, and 
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Hj = a f d^xe^^'n{x,t)(P{x,t) (B3) 

where 0(x, t) is a Hermitian operator made of any field or composite fields of the reservoir 
system that couples to the scalar field. For our derivation, we do not need specific details 
about Hr. It may well represent a strongly interacting system with respect to f2(x, t). 
Expressed in momentum space, eq. (Pl) is 



^3Ht 

Hi = -—Y,an{K,m-K,t). (B4) 

^ n 

Postulating the standard equal time CCR of eqs. and ([All|) , the equations of 

motion for the infiaton are 

0(k„, t) = i[Ht, 0(k„, t)] = e-3^*L37r^(k„, t) (B5) 



7r^(k„,t) = z[if,7r^(k„,t)] = -^[ul{t)<P{K.t) + + an{K.t)] (B6) 

where 



uj^it) = Ve-2^*k2 + m2. (B7) 
This implies the second order equation of motion 

AT/ 

0(k„, t) + 3i/0(k„, t) + (t)0(k„, t) + f + an{kn, t) = 0. (B8) 

d(f){-K,t) 

We derive the effective equation of motion for a given comoving mode of the infiaton. In 
the parametric regime eq. (|B1|), the time interval l/F ^ At ^ is sufficient to study 
the relaxational dynamics of the infiaton. Thus we can drop the second term in eq. ( P8|) 
above and assume that the comoving wavenumber k„ is constant at e~^^*k„ = k^, where 
k^ is the initial value in the time interval m of interest. Dropping the superscript m, the 
resulting second order equation of motion is 

6Vf 

0(k„, t) + ujI<P{K, t) + ; + «l](k„, t) = 0, (B9) 

where 

Un = + (BIO) 

This can be obtained from the first order equations in eqs. ( |B5| ) and ( |B6D by the replacement 

e-^^X _ e-^'^'-kl (Bll) 

in uJn{t) in eq. ( p6|) and by setting H = everywhere else. 

Doing this, the solution to eqs. (BSI) and (BBI) are respectively 
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0(k„,t) = 0o(kn,t) / di!VL{\<.n,i!) ^YnuJrSt - i!) 



- (B12) 



7r(k„,t) =7ro(k„,t) / (it'fi(k„, t') cosc<j„(t — t') 



- fdt' J^i^K. cosu;4t-t'), (B13) 
Jo d(p[-kn,t') 



where 0o and ttq are the solutions of eqs. ( P5| ) and ( |B^ ) with = and a = 0. Again 
for notational convenience we denote the physical momentum simply by k„. The associated 
comoving mode can be computed at time tm from eq. (pil|) . 

The inflaton (system) can be treated as a small perturbation on f2(k„,t), however full 
account must be taken of the heat reservoir's effect on the inflaton. Thus we will account 
for the lowest order effect of the system on the reservoir coordinate f2(k„,t), which gives 

fi(k„,t) = e-^^*fi(k„,0)e^^* 

icy \ 

- fdt'(i)iK',t-t')noi~K,t-t')noiK,t)\ (bi4) 

n' ^ ) 

where 

QoiK, t) = e-^^«*fi(k„, 0)e^^«* (B15) 

is the time development of f2 in the uncoupled reservoir. 

Substituting this into the second term on the right-hand-side of eq. (|B12| ), this term 
becomes 

— / dt'VL{kn,t') sinu;„(t - t') = — c/t'fio(k„, t') sincu„(t - t') + L(k„,t), (B16) 

UJn Jo UJn JO 

where 

L{kn,t) = —-Y, dt' sin Unit -t') dt" [<f){kn>,t")Qo{kn,t')Qo{-K',t") 
UJnL-^ ^ JO Jo 

- 0(k„,,t'-t")f^o(-k„',t'-t")^o(k„,t')]. (B17) 

We now take the statistical expectation value of eq. (p9|) with respect to the free states of 
the heat reservoir, which are assumed to obey a canonical distribution. The matrix elements 
of i7o(k„,t) between state < a\ and \b > of the uncoupled reservoir system are written as 

< a\noikn,t)\b >= {nf{kn,t)+znfikn,t))e-'^^^' (B18) 

where 
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(B19) 

with uj^ being the energy of the a-th state of the uncoupled reservoir. Hermeticity of 
imphes that Vtf' = and = — f^i". We assume that the diagonal matrix elements are 
zero 

n"^"" = 0. (B20) 

This is equivalent to assuming that there is no time independent force exerted by the heat 
reservoir on the system. The only nontrivial term for which the statistical expectation value 
needs to be evaluated is eq. (|B16|) . The expectation value of the first term on the right 
hand side of eq. ( |B16|) is zero due to eq. ( P20| ). However, we must retain this term in 
its operator form since we will also want products of it in the final rate equation when 
evaluating correlation functions. 

For L(kn,t) in eq. ( P16|) , upon statistical averaging we get 



N 

where 

2a 

ab 



^^(a|L(k„,t)|a)e--"/^ = L^(k„, t) + Lg(k„, t) (B21) 



L!ji{K,t) = ^—J2e-^/^Y. I dt' sin uj^{t-t') f dt" smujj^,{t' - t")<f){K' ,t")nl{a,b,K, 

UnlSL-^ ^ ^ Jo Jo 

(B22) 



^nNL^^ yjo Jo 

(B23) 



with 



and 



N = Y,e-</^, (B24) 

a 

nlia, b, K, k„o = (k„)^^^"(-k„o - r]f (k„)^]f (-k„o, (B25) 

nUa, b, k„, k„o = (k„)^]f (-k„o + Qf{K)Q',%-K'). (B26) 



The energy levels in the heat reservoir are assumed to be closely spaced so that we can 
replace the summation over the states by an integration as 

E - / P{u^^)du:^ / P(^f • (B27) 

We perform the above energy integration with the variables 
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= _ ^R^ 



(B28) 



so that 



POO POO roo roo rO roo p p 

j du^ j duj^ = j duo^j^Rdul + j duo^ j_^Rduo^ = j du^ j dCo^. (B29) 



Eqs. ( [B22|) and (|B23|) then become respectively 



00 



+ -f)pi^ 



~R 



dr 



cos(ci;:^^ + Un)^ sin(co':^ — Un) 



COs(co' 



~,R 



UJ- + UJn 



(k„,,t-r) (B30) 



1 



Define 



/ dr sin(co':^ + ujn) — sin(co':^ — cj^ 
Jo 2 



r 

'2 



(P{K,,t-T). (B31) 



roo 

Jo 



-uj' IT 



(B32) 



where I = C or I = D. We assume is a smooth function with respect to u. Identifying 
this term in eqs. ( [B30|) and ( [B31|) we observe that due to the oscillatory factor the main 
contribution to the uj^ integral is from uj^ = itu^. Evaluating the integration under this 
approximation, we obtain 



L^(k„, t) = Y^ r^„, /* dt'(f>{K,,t') cos Unit - t') 

I Jo 



(B33) 



where 



Lg(k„, t) = Y^ r^„, r dt'cj){K^,t') smcurrit - t') (B34) 



rL = ^^(1 - e--'^)B\un. k„, K') (B35) 
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with I = C^D. Substituting this into eq. (pi2|) we obtain the equation of motion 



(5 V 

0(k„,t) +u;2^(k„,t) + ^r^„,0(k„,,t) +^„^r^„,0(k„,,t) + = viK,t), (B36) 



where ?7(k„, t) = — af2o(kn, t). Hereafter we will assume that the uncoupled reservoir system 
is translationally invariant so that 

{{noiK,t)noiK', t')))T = 5n,-n'{{no{K,t)noi-K,t)))T. (B37) 

The statistical properties of ?7(k„,t) are 

{{r]{K,t)))T = (B38) 



((r^(k„, t)r/(k„,, t')))T = K,-n' + 7r<5(t - H) coth f ^ ) ) . (B39) 



The first property follows from eq. ( |B2CI| ). The second follows from eqs. ( [B25| ), ( |B26D 
( [B32D , and (^35]) and using 



OO np 

rfcje*-* = _ + vr5(t). (B40) 
t 



In the high temperature limit, eq. ( [B39| ) becomes 

((r/(k„,t)r/(-k„,t')))T^oo -2L3rJ„r<5(t-t'), (B41) 
which agrees with eq. 
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